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Inspirals of stellar-mass compact objects (COs) into ~ 10^ M© black holes are especially inter- 
esting sources of gravitational- waves for the planned Laser Interferometer Space Antenna (LISA). 
The orbits of these extreme-mass-ratio inspirals (EMRIs) are highly relativistic, displaying extreme 
versions of both perihelion precession and Lense-Thirring precession of the orbital plane. We investi- 
gate the question of whether the emitted waveforms can be used to strongly constrain the geometry 
of the central massive object, and in essence check that it corresponds to a Kerr black hole (BH). 
For a Kerr BH, all multipole moments of the spacetime have a simple, unique relation to M and 
S, the BH mass and spin; in particular, the spacetime's mass quadrupole moment Q is given by 
Q = —S'^/M. Here we treat Q as an additional parameter, independent of S and M, and ask how 
well observation can constrain its difference from the Kerr value. This was already estimated by 
Ryan, but for the simplified case of circular, equatorial orbits, and Ryan also neglected the signal 
modulations arising from the motion of the LISA satellites. We consider generic orbits and include 
the modulations due to the satellite motions. For this analysis, we use a family of approximate 
(basically post-Newtonian) waveforms, which represent the full parameter space of EMRI sources, 
and which exhibit the main qualitative features of true, general relativistic waveforms. We extend 
this parameter space to include (in an approximate manner) an arbitrary value of Q, and then 
construct the Fisher information matrix for the extended parameter space. By inverting the Fisher 
matrix, we estimate how accurately Q could be extracted from LISA observations of EMRIs. For 
1 year of coherent data from the inspiral of a lOM© black hole into rotating black holes of masses 
IQ^-^Mq, IO'^Mq, or IO'^ '^Mq, we find A{Q/M^) ~ 10"", 10"^ or 10"^ respectively (assuming 
total signal-to-noise ratio of 100, typical of the brightest detectable EMRIs). These results depend 
only weakly on the eccentricity of the inspiral orbit or the spin of the central object. 



I. INTRODUCTION 



Inspirals of stellar-mass compact objects (COs) captured by massive (~ IO^Mq) black holes (MBHs) in galactic 
nuclei are especially interesting gravitational- wave sources for the planned Laser Interferometer Space Antenna (LISA). 
Because the ratio of CO mass /i to MBH mass M is typically /i/M ^ 10~^, these events are generally referred to 
as extreme- mass-ratio inspirals (EMRIs). Current estimates are that LISA will detect of order ~ 10^ EMRIs per 
year P, [2, 3- Captures of ^ IOMq black holes (BHs) are expected to dominate the detection rate (although 
some captured white dwarfs and neutron stars will be detected as well), both because mass segregation in galactic 
nuclei tends to concentrate BHs very near the central MBH and because the inspiraling ^ lOM© black holes are 
intrinsically stronger sources than the white dwarfs or neutron stars, and so are detectable to a much greater distance: 
out to 2; ~ 1. 

The orbits of the inspiralling COs are highly relativistic, displaying extreme versions of both perihelion precession 
and Lense-Thirring precession of the orbital plane. Because of the tiny mass ratio, the inspiral is quite slow; individual 
inspiral waveforms will be observable for years, or equivalently for > 10^ cycles. While at any instant the EMRI 
waveforms will be buried in both EISA's instrumental noise and a gravitational-wave foreground from Galactic binaries, 
it will nevertheless be possible to dig out the EMRI signals using techniques based on matched filtering f3|, thanks 
to this large number of cycles. The matched-filtering signal-to-noise ratio (SNR), which scales like the square root of 
the number of observed cycles, will typically be in the range ~ 30 — 300 for EMRI detections. Because of this large 
SNR, and because the complex waveforms are quite sensitive to the source's physical parameters, it will be possible 
to extract these source parameters with high accuracy; in particular, it will be possible to infer the mass of the CO, 
and the mass and spin of the MBH, all to within fractional accuracy ^ 10~^ 

The famous "no hair" theorem of General Relativity essentially states that, almost immediately after its formation, 
a BH "settles down" to one of the Kerr solutions, and so its entire geometry is characterized by just two physical 
parameters: its mass M and its spin angular momentum S. In particular, all the multipole moments of a Kerr 
spacetime can be expressed in terms of M and S alone, according to the relation 



Mi+iSi = M{ia)\ 



(1) 
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where Mi and Si are the mass and mass-current multipole moments, respectively, and a = S/M is the standard Kerr 
spin parameter. For example, the mass quadrupolc moment Q = M2 of the pure Kerr geometry is given by 

Q = -S^/M. (2) 

It has been suggested by several authors (e.g., [1, 0, @|) that LISA detections of EMRIs could provide a highly 
accurate observational test of the "Kerr-ness" of the central massive object. (In what follows we will continue to 
refer to the central massive object as "MBH" for brevity, even though in our analysis we allow it to be a non-black 
hole object.) Indeed, Ryan Q showed that for inspiralling trajectories that are slightly eccentric and slightly non- 
equatorial, in principle all multipole moments of the spacetimes are redundantly encoded in the emitted gravitational 
waves, through the time-evolution of the orbit's three fundamental frequencies. (Basically, these are the fundamental 
frequencies associated with the 0, r and 9 motions; see ITo| for a proof that the waveform corresponding to a Kerr 
geodesic is a discrete sum of integer harmonics of these 3 fundamental frequencies.) 

Ryan [?! went on to estimate how accurately the central massive object's lower-order multipole moments could 
be measured by LISA. Because the CO's complex motion depends strongly on the background geometry, one might 
expect to be able to measure the lowest moments rather well. Basically, the long-term phase evolution of the orbit (and 
that of the emitted gravitational wave) can act as a microscope, probing the fine features of the geometry. Lacking a 
theory that accurately describes the orbital evolution in spacetimes with arbitrary multipole structure (in fact, such 
a theoretical framework is not yet at hand, even for the pure Kerr case), Ryan considered a simple waveform model 
based on lower-order post-Newtonian theory. For simplicity, he also restricted the EMRI parameter space to orbits 
that are circular and equatorial, and he neglected the modulations imposed on the observed signal by the spacecraft 
motions. Within these simplifications, he calculated the Fisher matrix on the space of EMRI waveforms, and inverted 
it to obtain (to lowest order in SNR^^) the variance-covariance matrix, which describes the distribution of statistical 
errors associated with a given measurement. Ryan found, for example, that for a capture of a lOM© black hole by 
a IQ^Mq central object, assuming 2 years of inspiral data with SNR of 100, LISA could measure Q/M^ to within 
~ 5 X 10~^, M to within a fractional error of - 3 X 10-4, S/NP to within - IQ-^, and ^ to within - 10~3. (Ryan 0| 
further showed how these accuracies get degraded as additional multipole moments, besides Q = M2, get treated as 
independent quantities in the signal model.) 

In essence, our paper improves on Ryan's 1997 analysis Q by considering the case of generic (i.e., non-circular, 
non-equatorial) orbits, and by employing a much more realistic treatment of EISA's instrumental response, including 
the Doppler modulation due to the constellation's center-of-mass motion and the additional modulations due to 
the rotational motion of the individual satellites with respect to that center-of-mass. Our analysis is based on our 
previously introduced family of "analytic kludge" (basically post-Newtonian) waveforms [H*], representing practically 
the full parameter space of capture sources (though neglecting the spin of the CO, for reasons described below). 
We extend this parameter space to include, in an approximate manner, an arbitrary deviation in the value of the 
central object's quadrupole moment, and we then construct the Fisher information matrix for the extended parameter 
space. By inverting the Fisher matrix, we estimate how accurately the central object's quadrupole moment Q can be 
determined, independently of its mass and spin. 

A priori, it would seem difficult to predict how our more realistic treatment should affect Ryan's results. On the 
one hand, our more realistic and complicated waveforms clearly encode more information. On the other hand, our 
waveform depends on all 15 physical parameters (including Q), while Ryan's more simplistic signal model was based 
on only 7 parameters j_29]; of course, one's ability to extract any given source parameter (such as Q) is diluted when 
more parameters are added to the source description. What we find, in the end, is that Q/M^ can be measured to 
^ 10"^ — IQ-^ (based on the last year of data and assuming SNR= 100). Thus EMRIs provide an even better test of 
"Kerr-ness" than Ryan estimated. 

Again, the idea of the Kerr-ness test explored here is to treat M, S, and Q as three independent quantities, and 
then see whether the LISA data confirms the Kerr relation: Q = —S'^/M. This test is similar in spirit to the well- 
known tests of general relativity made with the binary pulsar PSR 1913-1-16 (see, e.g., [HI). The binary's Keplerian 
parameters (like its period and eccentricity) and it post-Keplerian parameters (e.g., the precession rate uj, Einstein 
delay 7, and period derivative Pt) are all fit for independently. Within General Relativity these are redundant, since 
all parameters excluding Ph lead to a unique physical description of the binary. The test is then to see whether 
the measured Pb is indeed the value that GR predicts. Our EMRI test treats Q in an analogous way to Pi, for 
PSR 1913-1-16. Of course, testing that the spacetime's quadrupole moment coincides with the Kerr value is only 
one possible test of "Kerr-ness" , and not necessarily the best one. But absent a compelling alternative to the "null 
hypothesis" that the MBH is a Kerr BH, tests of this sort are perhaps the best one can do. 

We briefly mention some recent related work. In [l^], Collins and Hughes began a research program aimed at 
extending Ryan's work, beginning with the construction of exact spacetime solutions corresponding to "bumpy" 
black holes. Also recently, Glampedakis and Babak [l^ analyzed geodesic orbits in both Kerr and slightly non- 
Kerr spacetimes, and concluded that any effects due to non-Kerr multipole moments would be difficult to detect. 
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since they would be mimicked by small, compensating offsets in the geodesic's initial conditions. Glampedakis and 
Babak [12^ did not consider the effects of radiation reaction — i.e., they effectively considered only short stretches of 
actual waveforms — so the implications of their work for actual, inspiralling trajectories remained unclear. (Work to 
include radiation reaction in this model is now in progress (l4|). Finally, Kesden, Gair, and Kamionkowski jTS'] have 
investigated the question of whether EMRI waveforms can be used to distinguish between a central black hole and 
a massive boson star. They argue persuasively that one could easily distinguish between these two types of central 
object, since in the Kerr case the waveform "shuts off" shortly after the CO reaches the last stable orbit (LSO), while 
in the boson star case the inspiral continues until the CO comes to rest at the center of the boson star. 



II. "ANALYTIC KLUDGE" EMRI WAVEFORMS 



A. Overview 



Since highly accurate inspiral waveforms for EMRIs are not yet available (due to the difficulty of solving the 
radiation reaction problem for generic orbits in Kerr [16|), approximate, or "kludged", families of waveforms have 
been developed for use in EMRI investigations until more exact versions are available. One such family is the 
"numerical kludge" waveforms developed by Gair, Glampedakis and collaborators [13, [H, [3, and another is the 
"analytic kludge" waveforms developed by Barack and Cutler (hereinafter BC Q). Our "analytic kludge" waveforms 
are less accurate than the "numerical kludge" ones, but they are simpler to compute and more tractable for the 
purpose of accurately calculating partial derivatives (with respect to source parameters) numerically, as required for 
producing the Fisher matrix. The BC waveforms were used for just this purpose in Ref. [5(], where it was estimated 
that LISA could measure the masses of both bodies and the spin magnitude of the massive black hole to fractional 
accuracy ^ 10~^ — 10"'* (for SNR =100). In this paper we modify the analytic kludge evolution equations to include 
the influence of possible deviations in the quadrupole moment of the MBH. We now summarize this "Q-enhanced" 
version of our waveform family; we refer the reader to Ref. [5] for more details. 

In the "analytic kludge" approach, we approximate the CO-MBH system as being, at any instant, a Newtonian- 
orbit binary emitting a Peters-Matthews [23| (i.e., lowest order, quadrupolar) waveform. We then use post-Newtonian 
(PN) equations to secularly evolve the parameters of the orbit. In particular, we include orbital decay from radiation 
reaction, pericenter precession, and Lense-Thirring precession of the orbital plane. The motion of the LISA detector 
introduces additional modulations; for these we adopt the low- frequency approximation developed by Cutler [2]| . 
Cutler's treatment is not very accurate for frequencies / ^ 30 mHz (where the gravitational-wave period is longer 
than the round-trip travel time for light up and down the LISA arms), but it is adequate for our EMRI analysis, since 
most of the EMRI SNR accumulates in the range ^1 — 10 mHz. Since the main physical effects are all included, and 
their dependence on the source parameters are given correctly to (at least) the lowest non-trivial post-Newtonian order, 
we expect that parameter estimation accuracies for our "kludge" waveforms are a reasonable guide to the accuracies 
we can expect (i.e., are probably correct to within an order of magnitude), once the true, general relativistic waveforms 
are in hand. We test this expectation to some degree in Sec. III.C, where we simply remove (by hand) some P^N 
terms from our kludged evolution equations, and find that this modifies the deduced measurement error in Q /A/P by 
a factor of less than 3. 



B. Parameter Space 

Assuming the MBH is a Kerr back hole, the two-body system is completely described by 17 parameters. However, 
the spin of the CO is at most marginally relevant so, as in BC, we shall neglect it in our analysis, [s^ That leaves 
14 parameters, assuming a massive Kerr BH. We now include one additional parameter, describing the quadrupole 
moment of the geometry. More precisely, we choose as our additional parameter the dimensionless quantity 

Q = Q/Af3. 

Including Q, we therefore have a 15-dimensional parameter space. For completeness, we list the full parameter set 
here. This list is identical to the one in BC, except for the inclusion of Q (and a re-definition of to — see below). 

A"^ = (A",...,A") = 

to{xl niHz), In^, InM, S, clso, 7o, *o, Ms = cos^s, 413, cos A, ao,^J.K = cos6k, (f>K, ln{n/D),Q . (3) 

Here, to specifies the instant of time where the CO reaches the LSO, and the inspiral transits to a plunge, fi and 
M are the masses of the CO and MBH, respectively, and S = S/M'^, where S is the magnitude of the MBH's spin 
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angular momentum. The parameters Clso, 7o, and $o describe, respectively, the eccentricity, the direction of the 
pericenter within the orbital plane, and the mean anomaly — all at time tg (time of plunge). More specifically, we 
take $0 to be the mean anomaly with respect to pericenter passage, and 70 to be the angle (in the plane of the orbit) 
from L X S to pericenter (where L is a unit vector in the direction of the orbital angular momentum). The parameter 
ao = a{t = to) describes the direction of L around S at time to- [See Eq. (18) of BC for the precise definition of 
a{t).] The angles {9s,4>s) describe the direction to the source, in echptic-based coordinates; {6k, (px) represent the 
direction S of the MBH's spin (approximated as constant) in ecliptic-based coordinates; and A is the angle between L 
and S (also approximated as constant). Finally, D is the distance to the source, and, again, Q = Q/M^ is spacetime's 
dimcnsionless quadrupole moment. The various parameters and their meaning are summarized in Table [D Fig. 1 of 
BC illustrates the various angles involved in our parameterization. 

Note that for simplicity we are treating the background spacetime as Minkowski space, not Robertson- Walker. 
To correct this, for a source at redshift z, requires only the simple translation: M M{1 + z), ^ ^ ^(1 + z), 
S S{1 + z)^, D Dl, where Dj^ is the "luminosity distance" |23 |. 



A" 


to(x ImHz) 


time when CO reaches LSO 




\nfi 


(In of) CO's mass 


A^ 


InM 


(In of) MBH's mass 


A^' 


S = S/M^ 


dimcnsionless magnitude of (specific) spin angular momentum of MBH 


A* 


fiLSO 


eccentricity at % (i.e., final eccentricity) 


A^ 


70 


7(io), where 7(t) is the angle (in orbital plane) between L x S and pericenter 


A« 


$0 


$(fo), where $(f) is the mean anomaly 


A^ 


/is = COS Os 


(cosine of) the source direction's polar angle 


A« 


(ps 


azimuthal direction to source 


A^ 


COS A 


L ■ S{= const) 


Aio 


QO 


a{to), where a{t) is the azimuthal direction of L (in the orbital plane) 


A" 


/iK = COS Ok 


(cosine of) the polar angle of MBH's spin 


A^^ 


4>K 


azimuthal direction of MBH's spin 


A" 


Hp/d) 


(In of) CO's mass divided by distance to source 




Q = QIM"" 


dimensionless quadrupole moment of MBH 



TABLE I: Summary of physical parameters and their meaning. For further details see Fig. 1 of BC and the description in the text. 



C. Orbital evolution with arbitrary quadrupole moment 

The inspiralling orbit is determined by the evolution equations for ^{t),v{t),^{t),e{t),a{t). (The 2 masses, the 
A/IBH spin vector 5*, the angle A and of course the source's sky position are all taken to be constants.) The post- 
Newtonian evolution equations we adopt are: 



2ttv, (4) 



^\ = —[^Ji/M^){2nMvf^'\l - e2)-9/2{ [1 + (73/24)e2 + [il /m)e'] (1 - e^) 
at J IOtt 

-|-(27rAfiy)2/3 [(1273/336) - (2561/224)6^ - (3885/128)e'* - (13147/5376)e^] 
-(27rAfiy)S' cos A(l -e2)-i/2[(73/;L2)-f (1211/24)6^ + (3143/96)6-^-)- (65/64)e^] 
-{2T:Mvfl^Q{l - e^)-! [(33/16) -I- (359/32)6^ - (527/96) sin^ A] }, (5) 

^ 6^i.(2^i.A/)2/3(l-e2)-i l + i(27ri.M)2/3(l_e2)-i(26-15e2) 
-127riycosAS'(27rMi/)(l - e^)-^/'^ 

-^7riyQ(27rMi/)'*/3(l-e2)^^(5cosA- 1), (6) 



5 



de 



dt / 15 



- — {^l/M'^){l ~ e2)-7/2(27rAf,.)8/3[(304 + 121e2)(l - e^)(l + l2{2TiMvf/^) 



-^(27rMi/)2/3((8)(16705) + (12)(9082)e2 - 25211e'')] 

+e(^/M2)S'cos A (27rAfiy)"/3(i _ ^^yA [(1364/5) + (5032/15)e2 + (263/10)e''] , (7) 
^ 4^z.5(27rMz.)(l-e2)-3/2 

+iTivQ{2'KMvf'^{l - e^y'^cosX. (8) 

These are the same evolution equations as used in BC, except that the lowest-order terms cx Q have been added to the 
right-hand side of Eqs. (O, ([6]), and ([8]). For the two precession rates, dj/dt and da/dt, these extra terms represent 
essentially Newtonian effects: the Newtonian precession rates for orbits around oblate bodies. The term cx Q in our 
Eq. ^ is taken from Eq. (3) of Lai et al. [2^, using the identifications X ^ 7, /i — /a — > Q and 6 \, and the 
approximation k, 9c for ji/M <^ 1. The term oc Q in Eq. ([5]) is taken from Eq. (4) of Lai et al. To obtain the 
term cx Q in Eq. ([5]) for we start with the Newtonian relation 

^ {-2E/t,f/\ (9) 



27rM 

This relation is not modified by Q, through first order [cf. Eq. (2.26) of Gergely et al. [l^l or Eqs. (A.11)-(A.12) of 
Ashby [25|]]. Therefore we can write, through first order in Q and to lowest non-trivial order in a PN expansion: 

yQ = = --y{2T:Myr^'^EQ . (10) 

CLrj H 

where the subscript 'Q' here means "the piece linear in Q". For Eq, we adopt the approximate formula given in 
Eq.(44) of Gair and Glampedakis [l^. [As those authors note, their Eq. (44) is unfortunately "missing" a piece 
proportional to sin'^ A, since the coefficient of that term has not yet been derived.] Using the identifications that 
their corresponds to our —Q and their t to our A, and then plugging Eg into Eq. pI7|) . we obtain the term cx Q in 
Eq. ©. 

For consistency, our Eq. ([7]) should also contain a term cx Q; however it would be nontrivial to derive this term 
from the existing literature, and its influence on astrophysically realistic waveforms is probably substantially less than 
the terms cx Q that we have included, since ^ 0.1 for typical inspirals. It is a feature of our kludged waveforms 
that we attempt to include the most important terms but do not worry about consistently including all terms out to 
some fixed order; in this spirit, we simply neglect the Q-correction to e. 

To construct the orbit, we start by choosing some final eccentricity Clso- In practice we consider moderate final 
eccentricities: those in the range < clso < 0.3. We then set the final radial frequency as 

t-Lso = 0.0415 X {2'kM)-\ (11) 

which reproduces the correct gravitational- wave frequency at the LSO through Eq. ([H)) (explained further below). 
This LSO value of the frequency applies to circular orbits in Schwarzschild (5 = 6 = 0), but for simplicity we shall use 
it here for any of the moderate eccentricities considered, and for any spin. We then integrate the evolution equations, 
Eqs. (4)-(8), backwards in time, starting from the LSO. Given the solution of these ordinary differential equations, it 
is straightforward to construct both polarizations of the emitted waveform using Eqs. (7)-(10) in BC: 

h+,At) = Yl [<' " ^i'^*" (^) + ' COS*" (^)] ' (12) 

n 

where the subscript "n" labels the different harmonics of the orbital frequency, and the subscripts "s" and "c" just 
refer to the two quadratures, sine and cosine. In practice, we find it sufficient, given the range of LSO eccentricities 
considered in this work, to include in the above sum only the first 20 harmonics. Explicit formulae for the and 
A^'^ are trivially obtained from Eqs. (7), (8), and (10) of BC, and the phase of the n*'' harmonic, <I>„(t), is given by 

$n(i) = n$o + 27r /" fn{t')dt' + 27r/„(t)n • R{t) . (13) 

Jto 
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Here /n(i') is the instantaneous frequency of the n*^ harmonic, and the second term in Eq. (fT3|) is the Doppler delay, 
where n is the direction to the source and R{t) is the position of LISA's center with respect to the Solar System 
Barycenter [3l|. For our "kludge" waveforms, we approximate /„(i') by simply 

Mt') = n,y{t')+^{t')/n, (14) 

the motivation for which is explained in the last paragraph of Section IV of BC. 

To explain our choice for t'Lso [Eq- (|lip above], recall that, for circular orbits around a Schwarzschild black hole, the 
(azimuthal) orbital frequency at the LSO is precisely v^so = 6~'^/^(27rM)~^, with a corresponding gravitational wave 
frequency of /gw = 2 x 6~'^/^(27rAf)~-'^ (since the radiation in this case is dominated by the n = 2 harmonic). To be 
able to reproduce this value of /gw at the LSO through Eq. ([M]) . we must require 2fLSO+7LSo/''' = 2x 6^^/^(27rAf )~^, 
where 7lso is obtained from Eq. ([6]) by replacing v — > 14,30 and neglecting the terms proportional to e^, S, and Q. 
Solving this condition for ^lso gives Eq. pT|) above. 

In the frequency range of interest, we can regard LISA's output as equivalent to two independent data streams, si{t) 
and sii(t). We idealize each of these as being the sum of (stationary, Gaussian) noise niji{t) and a gravitational- wave 
signal, hiji{t): 

si{t) =ni{t) + hi{t), .3iiit)^nii{t) + hii{t). (15) 

The gravitational- wave signals in the two channels, hi{t) and hii{t), are related to the two incoming polarizations 
h^^x{t) by time-varying transfer functions, F^'^{t) and F^j'^{t). Thus we have 

hiit) - y^[F+{t)h+it)+F,-it)hxit)] , (16) 

and similarly for hii{t). The equations used to construct these transfer functions are given explicitly in BC. (The 
factor V3/2 is just the sine of 60°, the angle between the LISA arms; it could have been included in the definition 
of the transfer functions, but it was left as a separate factor so that the F's would coincide with LIGO transfer 
functions.) 

We emphasize that in our "analytic kludge" scheme, the value of the quadrupolc moment affects the inspiral 
trajectory (i.e., it modifies the sequence of Keplerian ellipses that the orbit osculates though), but Q does not enter 
into the algorithm for constructing the gravitational waveform from that trajectory; i.e., we continue to use the 
quadrupole formula to construct the waveform from the inspiral trajectory. Since the quadrupole formula does a 
surprisingly good job of reproducing waveforms for Kerr [l^, we expect this also to be the case for spacetimes that 
are slightly off-Kerr. 

Finally, we briefly review the signal analysis formalism necessary for calculating AQ (the measurement accuracy 
in Q) in order to make clear our approximations; we refer to BC for more details. The (sky-averaged) noise spectral 
density Sh{f) for each of the two channels, I and II, determines a natural inner product on the space of LISA data 
streams, as follows. Let p(i) = [p/(i),p//(t)] and q(t) = [qi{t)^qij{t)] be two different LISA data sets. Their inner 
product is [26l|: 



/"OO r q 

(p|q)^2^ / WM)Uf)+PMWM)\l 



df . (17) 



[The factor ^ is due to the "sky-averaging" convention adopted here in the definition of Sh{f).] This inner product 
has the property that the matched-filtering SNR for any imbedded gravitational waveform h(<) = [hi{t),hjj{t)\ is 
just its norm: 

SNR[h] = (hjh)'/' . (18) 
The Fisher information matrix Tab is defined by 



Fab 



dh 



(19) 



where a, 5, . . . are parameter-space indices. To lowest order in an expansion in SNR ^, the variance-covariance matrix 
for the errors is just the inverse of the Fisher matrix: 

{SrSX'') = (r-1)"^ [1 + ©(SNR-i)] , (20) 
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where "(• • •)" means "expectation value". Hence the measurement error in the various parameters is given by 

AA'^ = {SX'SV)^''^ w [{T-^yY/'^ . (21) 

The actual inner product, Eq. ()17|) . is formulated in the frequency domain. For a white noise [i.e., Sh{f) ~ constant], 
the inner product is equivalent to 25*^^ ^^.^j/^ Pa{t)qa{t)dt, by Parseval's theorem. Motivated by this formula, 
we shall adopt the following approximate version of the inner product in calculating the Fisher matrix: First, we 
define the "noise-weighted" waveform 

h^{t) = J2hcAt)/Sl/'{fn{t)), (22) 

n 

where (again) we take fn{t) — ^{t) -\-j{t)/7T. Then we approximate the covariance matrix, Eq. (jl9p . as 

rafc = 2(|-)V/ daK{t)dM)dt. (23) 

That is, we simply re-weight each harmonic by the square root of the inverse spectral density of the noise, and 
thereafter treat the noise as if it were white. 

Finally, we must specify the (sky-averaged) noise spectral density Sh{f) that we use in constructing the inner 
product. For this, we use the noise model given in Eq. (30) of |27|]. This model treats the noise as the sum of 
instrumental noise and confusion noise from galactic and extragalactic white dwarf binaries. All low- frequency binaries 
contribute additively to the confusion noise, but the effect of the high-frequency galactic binaries (which can all be 
detected and fitted out of the data) is to effectively reduce the bandwidth for observing other sources (like EMRIs) by 
some factor, or equivalently to increase the effective noise density Sh{f ) by that same factor; see [27j for more details. 



III. RESULTS 

In this section we present our results for AQ, the accuracy with which LISA can determine the MBH's (dimension- 
less) quadrupole moment Q = Q/M^ , for a range of physical parameters. This paper does not aim at a full, Monte 
Carlo-type exploration of AQ over the entire 15-dimensional parameter space. Rather, we regard this paper as an 
initial exploration, aimed at determining the order of magnitude of AQ for a few "typical" cases. Based on this goal, 
our survey strategy was as follows. We first attempted to find values of the angular parameters {9s,4>S,()k,4'k, ^) 
that led to median results for AQ. We then fixed these angles and explored the rest of parameter space. We always 
evaluated the Fisher matrix at the value of Q corresponding to pure Kerr geometry, i.e., Q = — (S'/M^)^; thus, in 
effect, we explored how accurately small deviations from the Kerr value could be measured. In Sec. IIII A1 we describe 
how we obtained "representative" values for the angular parameters. In Sec. IIII B] we present our results for AQ, for 
a variety of values of M, fi, S, and blso- Finally in Sec. IIII CI we present checks on the accuracy and robustness of 
our numerical results. 



A. Choice of Representative Angles 

Here we describe our method for choosing one set of representative angles. As we show below, AQ depends strongly 
on M, and more weakly on ^, clsoj and S. We therefore fixed /i = IOMq, clso = 0.15, and S = 0.5, and for each 
of the MBH masses M = IO^Mq, and IO^-^Mq explored all 108 different values of {es,(t>s,dK,(t>K,X) 

from within 6*5 G (^tt, ^tt, |7r) , (j}s £ (O, |7r, |7r), 9k e (^^r, ^tt, |7r) , (f>K £ (O, ^tt, tt), and A € ( j^, ^tt, |7r) . We 
calculated AQ for each of these points and found the median values of AQ for each of M = 10^'^, 10^, 1O^'^M0. The 
results for M — IO^Mq are histogrammed in Fig. [1] The following choice of angles turns out to yield a value of AQ 
close to its median: 

[ (f,f ,f,0,f), for Af = IO^-SMq, 
{esAs,dK,^K,X) = I (^,3^,f,0,f) , forM^lO^MQ, (24) 
[ (|,f ,f ,7r,f), for Af = IO^-SMq. 

We adopted these as our representative angles in surveying the (^, S, blso) space (for each of the different MBH values 
explored) . 
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mean[log^i,A (Q/M^)] = -2.964 (SVD) 
= -2.932 (Gauss) 

STD[log,^A (Q/M^)] = 0.130 (SVD) 
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FIG. 1: Distribution of AQ for a selection of (6*3, (jjs, Ok, 4>k, A) values, and fixed values of (/^, M, eLSo) = (IOM0, 10® Mq, 0.15). 
The mean, standard deviation (STD), and median of the distribution are calculated twice, for comparison, based on data from 
two different matrix-inversion methods: singular-value decomposition (SVD) and Gauss- Jordan elimination (see discussion in 
the text). We find that the choice [Os, 4's,dK, 4>k, A) — (^,^,-|,0, -1) yields AQ very close to the median value. The median 
values for the masses M = lO'^-^ and IO'^ '^Mq were obtained in a similar manner, and are given in Eq. (|24|l . The standard 
deviations in logj^Q AQ were found to be ~ 0.13 and ~ 0.17 for M = lO^ '^ and 10'^ '^ Mq, respectively. 



B. LISA'S measurement accuracy for Q 

In this subsection we present our results for AQ, obtained by calculating the Fisher matrix and using the leading- 
order estimate (AQ)^ = (F"^)'^'^. All results are based on analyzing exactly one year of data, which we take to be 
the last year of inspiral, and all are normalized to a total SNR of 100. Our reason for normalizing to SNR= 100 
is that the best tests of "Kerr-ness" will naturally come from the strongest (generally closest) EMRI sources. The 
detection threshold for EMRIs, as limited by computational power, is estimated to be SNR~ 30 — 35 Assuming 
> 100 EMRIs are detected above this threshold, an estimated > 4 EMRIs wiU be detected with SNR> 100. Thus 
SNR= 100 corresponds to a relatively strong source, but we do expect to detect a few sources that strong. Of course, 
AQ roughly scales like SNR~^, so even for marginally detectable sources with SNR^ 30 — 35, AQ will only be a factor 
~ 3 larger. 

Fig. [2] presents our results for AQ for two fiducial CO masses (lOAf0 BH and O.GMq white dwarf), three fiducial 
MBH masses (10^-^, 10^ °, and lO^-^M©), three fiducial values of clso (0, 0.15, and 0.3), and five fiducial values of 
S (0, 0.25, 0.5, 0.75, and 1). Again, all results in these figures are for the sets of "representative" angles specified 
in Eq. ([24]) . We have not included results for blso) = (lOM©, 1O^'^M0, 0.3) [and also missed out the point 

{fi, m, Chso, S) — (lOAf0, lO^'^M©, 0.3, 0)] as for these parameter values the initial eccentricity becomes too high 
(e ^ 0.7), and our evolution model fails to yield reliable results. (The failure of our kludge model at very high 
eccentricities is discussed in Appendix B of EC.) 

The following general conclusions are evident from the above Figures. For captured IOAIq BHs it should be possible 
to extract Q to within an accuracy AQ ~ 10~'*-10^^. The value of AQ varies strongly with M, being ^ 100 times 
larger for M = than for M = lO^-^M©. As ^JL decreases from to 10 to O.eM© (with SNR assumed fixed), 

AQ increases by a factor ^ 10. These patterns can be understood qualitatively as follows. At fixed SNR, one expects 
AQ, A(lnAf), A(ln/i), and A5 to all decrease with an increasing number of cycles, and the latter scales roughly as 
'^Lso (for a fixed observation time), or equivalently as . One also expects AQ, A(lnAf), A(ln/i), and AS' to all 
decrease with an increasing ratio of /2,LSo//2,yr, where /2,yr is the frequency of the n = 2 harmonic one year prior 
to plunge, and /2,lso is that harmonic's frequency just prior to plunge. That is because a larger ratio of final to 
initial frequencies provides more "leverage" for separating out post-Newtonian effects that enter at different orders of 
the FN expansion parameter, (7rM/2)^/^. We can obtain an approximate expression for /2,LSo//2,yr by simplifying 
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FIG. 2: Results for AQ — the measurement accuracy in Q = QjM'^ — for a selection of values of M, blso, and 5* = S/M'^ . 
Left and right panels display AQ for IOA/0 and 0.6A'/q COs, respectively. Both plots are for one year of data (the last year 
of inspiral), and are normalized to SNR= 100. [In the left panel we have discarded a few of the data points, corresponding to 
low MBH mass with high LSO eccentricity: For these points the initial eccentricity exceeds 0.7, and our PN evolution model 
cannot be trusted (see discussion in the text).] 



to the case of e = orbits in Schwarzschild, so /2,lso — (6^/^7rM) and by approximating df2/dt using the lowest 
(nontrivial) order PN expression: 



5 

which yields 



Qft /-lyr /•/2.LS0 

n^'/^liM^'^ dt^ r^^'^'df, (25) 



Jf2.yr 



/2,yr ; ■ lOi\/0 V M 



Therefore, for the /i — lOM© case, as M increases from 10 to 10 A/q, /2,LSo//2,yr decreases from 3.4 to 1.09. 
Similarly, for M = IQ^Mq, /2,LSo//2,yr decreases from 1.6 to 1.05 as decreases from 10 to O.GAf©. 

When one adds Q to the list of parameters that must be extracted, the accuracy with which all other parameters 
can be extracted must decrease. The magnitude of this effect, in our case, is illustrated in Table |TT1 where we list 
A(ln/i), A(lnAf), and AS* for our three fiducial M values and /i — IQMq. The bold-faced entries are the error 
magnitudes for the 15-dimensional parameter set including Q, while the plain-faced entries are for the 14-dimensional 
set that excludes Q. As expected, the errors when Q is included are always larger, by factors that range from 1.13 
to 54. Including Q has the greatest impact on AS*. This is also to be expected, since the terms oc Q and those oc S 
in Eqs. (O, (O, and ([5]) have a very similar scaling with frequency: the former are multiplied by {21: Mv)'^/'^ and the 
latter by {2'kMv); thus one expects errors in the two quantities to be strongly correlated, and so also for Q to have 
the strongest effect on S. The effect of the extra parameter Q on the size of A5 is largest for high A/, since the 
smaller the ratio /2,LSo//2,yr, the more degenerate these parameters become. 



C. Inversion checks and test of robustness 



The Fisher matrices we compute typically have a very large condition number (the ratio of the largest to smallest 
eigenvalues). For the numerical computation of a matrix inverse to be reliable, using double pre cision arithmetic, it 
is generally sufficient (but not necessary!) that the matrix's condition number be ^ 10^^ |2^. |32| Some of the Fisher 
matrices we calculate have condition numbers as large as ~ 10^^. This is cause for some concern when we numerically 
invert the Fisher matrices to find AQ; specifically, one might worry that small errors in the Fisher matrix elements 
(e.g., due to round off) could lead to large errors in the matrix inverse. Therefore we have performed several checks 
that our results are reliable. 
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M 


A(ln^i) 


A(lnAf) 


AS 




2.9 X 10" 

1.8 X 10~^ 


6.9 X 10 

4.7 X 10-^ 


8.3 X 10-^ 
1.7 X 10"'' 




5.2 X 10"^ 

2.5 X lO""' 


1.6 X 10" 

1.1 X 10"^ 


1.6 X 10 " 

3.6 X lO"'* 




8.8 X 10"'' 

5.9 X 10-" 


2.6 X 10" 

2.3 X 10"® 


2.6 X 10" 

4.8 X 10"^ 



TABLE II: Bold-face data in this table show the measurement accuracy of In fi, In Af, and S when Q is included as a search parameter. 
(Here We have taken ^ = WMq, blso = 0.15, S = M^, and a few values of M, which corresponds to three of the data points in Fig. 2 
above.) For comparison, plain-face data show the measurement accuracy of these parameters for the same points in parameter space, but 
without searching over Q, i.e., assuming a Kerr black hole: Q = —{S/M'^)"^. 



First, we have used several different matrix-inversion routines, and checked that the results from different routines 
are always fairly close together. This is illustrated in Fig. [31 where for the 108 matrix inversions used to generate Fig. 
[3 we plot the values of AQ derived from two different methods: (i) Gauss-Jordan elimination and (ii) singular-value 
decomposition [28^. The results from the two different inversion routines agree to within less than ^ 20% in all cases. 
Additional evidence that our numerical inverses are numerically stable come from the plots in Fig. [21 which show 
that AQ depends in a rather smooth and consistent way on changes in S and clso; i-e, the matrix inverse is clearly 
not wildly sensitive to changes in source parameters. Thirdly, if we recondition our matrices by essentially re-scaling 
all the variables, so that diagonal elements of the new Fisher matrix are all unity, then a matrix with an original 
condition number of ~ 10^^ typically has a new condition number as low as ~ 10^^. 

It should be pointed out that the error in some of the orbital phases (those associated with the smallest singular 
values) are strongly dependent on the inversion method, and our code does not produce reliable answers for these 
(apart from telling us that these parameters are strongly degenerate). Fortunately, the parameter Q is among those 
for which the inversion gives robust answers. 
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FIG. 3: Comparison of the values of AQ obtained using two different matrix inversion methods. The 108 sample points 
are the same ones displayed in Fig. [1] i.e., they correspond to 108 different values for {6s,4>s,0k,4'k,X), and fixed values 
of (/i, Af, blso, Q)- Circles are results using singular- value decomposition, while pluses correspond to inversion of the same 
matrices using Gauss- Jordan elimination. Generally the two results lie close together, with deviations smaller than ~ 20% in 
all cases. This is sufficient agreement for our purpose. 

A related worry is that our results for AQ might depend strongly on the exact signal model, so that AQ might 
be very different for our "analytic kludge" waveforms than for true, general relativistic waveforms. We now present 
evidence that, to the contrary, our general conclusions seem rather robust. Roughly speaking, our "analytic kludge" 
waveforms differ from the true ones due mainly to all the higher-order FN terms that are missing from the former. 
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(Our treatment of the LISA response is also approximate, but our low-frequency approximation here seems much 
less likely to significantly bias the final results.) We can get some sense of the importance of these missing terms 
for the final conclusion, by comparing our results with a "lower-order" or "dumbed-down" version of the analytic 
kludge. Specifically, the dumbed-down version just neglects all the correction terms proportional to (27rMz/)^/'^ in 
Eqs. ©-([S]) . That is, we simply drop the second line in Eq.(5), the second term in square brackets in the first line 
of Eq.(6), and the second line of Eq.(7). The remaining equations still have the property that the major qualitative 
features of the orbits (perihelion precession, Lense-Thirring precession, and inspiral) are all present, and that at least 
the lowest-order PN term corresponding to each effect is included. 

Our results for AQ, using this "dumbed-down" version of the analytic kludge evolution equations, are demonstrated 
in Fig. [3]for the case fi = IOMq. The values of AQ remain the same to within a factor ^ 2 (compare with the left 
panel of Fig. ^ . We take these results as indication that our estimates for AQ are accurate to better than an order 
magnitude — perhaps even to within a factor ~ 3. 
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FIG. 4: AQ for BH inspirals, for a lower-order version of our "kludge" evolution equations (given explicitly in the text). 
Comparison of this Figure with the left panel of Fig. [2] provides a check on the robustness of our results. 



IV. CONCLUSIONS AND DISCUSSION 



Our main conclusion is that, for reasonably strong EMRI sources (SNR^ 100), it should be possible to extract 
the quadrupole moment Q of the central massive object, independently of its mass and spin, typically to within 
AQ ~ 10~'*-10~^. Thus for each detected EMRI, LISA will have the opportunity to perform a rather precise and 
nontrivial check of the "Kerr-ness" of the central massive object. For fixed SNR, smaller values of M will provide 
much more accurate measurements of the MBH's quadrupole moment. While the particular source orientation and 
sky direction can cause AQ to vary by up to a factor 4, AQ seems to depend rather weakly on clso and on the value 
of the spin S. 

An obvious caveat is that the values of AQ derived here represent the magnitude of the error due to noise alone; any 
errors in the construction of theoretical waveforms for given parameter values would lead to additional, systematic 
errors in Q. Another caveat is that even if the best-fit values of Q were substantially off-Kerr in a few sources, this 
could just mean that those sources were not sufficiently "clean" two-body systems, e.g., because of the influence of 
an accretion disk or because of interactions with a third body. 

Finally, it might be interesting to extend our analysis to consider independent estimates of the first N multipole 
moments (for = 3, 4, 5, . . .). Indeed this is what Ryan [7,] did, within his highly simplified model. However we have 
not yet attempted this. 



12 



V. ACKNOWLEDGMENTS 

We thank for Kip Thorne for encouraging us to work on this problem. We thank Neil Cornish for pointing out a 
typographical error in Eq. (33) of BC. L.B.'s work was supported by PPARC through grant number PP/DOOlllO/1. 
L.B. also gratefully acknowledges financial support from the Nuffield Foundation, and thanks Kostas Glampedakis for 
helpful discussions. C.C.'s work was carried out at the Jet Propulsion Laboratory, California Institute of Technology, 
and was sponsored by the National Aeronautics and Space Administration. 



[1] C. Hopman and T. Alexander, Astrophys. J 645, L133 (2006). 

[2] M. C. Miller, M. Freitag, D. P. Hamilton and V. M. Lauburg, Astrophys. J 631, L117 (2005). 
[3] S. Sigurdsson, Class. Quant. Grav. 20, S45 (2003). 

[4] J. R. Gair, L. Barack, T. Creighton, C. Cutler, S. L. Larson, E. S. Phinney and M. Vallisneri, Class. Quant. Grav. 21, 
S1595 (2004). 

[5] L. Barack and C. Cutler, Phys. Rev. D 69 082005 (2004). 

[6] F. D. Ryan, Phys. Rev. D 52, 5707 (1995). 

[7] F. D. Ryan, Phys. Rev. D 56, 1845 (1997). 

[8] S. A. Hughes, Class. Quant. Grav. 18 4067 (2001). 

[9] W. Schmidt, Class. Quant. Grav. 19, 2743 (2002). 
[10] S. Drasco and S. A. Hughes, Phys. Rev. D 73 024027 (2006). 
[11] D. Lorimer, Living Reviews in Relativity (2005); astro-ph/0511258' 
[12] N. A Collins and S. A. Hughes, Phys. Rev. D69, 124022 (2004). 
[13] K. Glampedakis and S. Babak, Class. Quant. Grav. 23, 4167 (2006). 
[14] K. Glampedakis, private communication. 

[15] M. Kesden, J. Gair and M. Kamionkowski, Phys. Rev. DTI, 044015 (2005). 
[16] E. Poisson, Living Rev. Relativity 7, 6 (2004). 

[17] K. Glampedakis and D. Kennefick, Phys. Rev. D66, 044002 (2002). 

[18] J. R. Gair and K. Glampedakis, Phys. Rev. D 73, 064037 (2006). 

[19] S. Babak, H. Fang, J. R. Gair, K. Glampedakis, and S. A. Hughes, gr-qc/0607007 

[20] P. C. Peters and J. Mathews, Phys. Rev. 131, 435 (1963); P. C. Peters, Phys. Rev. 136, B1224 (1964). 
[21] C. Cutler, Phys. Rev. D. 57, 7089 (1998). 
[22] D. Markovic, Phys. Rev. D 48 4738 (1993). 

[23] D. Lai, L. Bildsten and V. M. Kaspi, Astrophys. J 452, 819 (1995). 
[24] L. A. Gergely and Z. Keresztes, Phys. Rev. D 67, 024020 (2003). 

[25] N. Ashby, Proceedings of 33rd Annual Precise Time and Time Interval Systems and Applications Meeting (Long Beach, 

CA), 2002; ]http://tycho.usno.navy.mil/ ptti/ptti2001/pape r51.pdf| 
[26] C. Cutler, and E. E. Flanagan, Phys. Rev. D 49, 2658 (1994)7^ 
[27] L. Barack and C. Cutler, Phys. Rev. D 70 122002 (2004). 

[28] W. H. Press, S. A. Teukolsky, W. T. Vetterhng, and B. P. Flannery, Numerical Recipes (Cambridge University Press, 
Cambridge, 1992). 

[29] Our additional 8 parameters are: 2 sky-position angles, 2 angles describing the MBH's spin direction, 2 angles describing 
the initial orientation of the orbital angular momentum, 1 angle describing the direction to pericenter, and the final value 
of the eccentricity, clso- 

[30] More specifically, it was shown in Appendix C of 5] that effects related to the CO's spin can modify the waveform phase 
by at most a few radians out of ^ l(f , for a maximally spinning CO. This is far smaller than the sought-for influence of an 
excess quadrupole moment for the central object: A quadrupole moment twice its Kerr value would modify the year- long 
inspiral waveform by ~ 10^ radians. 

[31] Eq. (33) of BC incorrectly gives the Doppler piece of "l>„(t) as 2-Knv{t)n ■ R{t). However the code that generated the results 
in BC used the same approximation as given in this paper: 2nf„{t)h ■ R{t). 

[32] As an example of a matrix that fails this condition but that is easy to invert numerically, consider the 2x2 matrix with 
diagonal elements 1.0 and 10""^", with zeroes for the off-diagonal elements. 



